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Abstract The planar quantum dynamics of spin-1/2 neu¬ 
tral particle interacting with electrical fields is considered. A 
set of first order differential equations are obtained directly 
from the planar Dirac equation with nonminimum coupling. 
New solutions of this system, in particular, for the Aharonov- 
Casher effect, are found and discussed in detail. Pauli equa¬ 
tion is also obtained by studying the motion of the parti¬ 
cle when it describes a circular path of constant radius. We 
also analyze the planar dynamics in the full space, including 
the r = 0 region. The self-adjoint extension method is used 
to obtain the energy levels and wave functions of the par¬ 
ticle for two particular values for the self-adjoint extension 
parameter. The energy levels obtained are analogous to the 
Landau levels and explicitly depend on the spin projection 
parameter. 

1 Introduction 

Topological effects in quantum mechanics has been one of 
the most studied problems of planar dynamics in recent years. 
These phenomena present no classical counterparts and are 
associated with physical systems defined on a multiply con¬ 
nected space-time [1]. Many of the recent interest in this 
matter is a consequence of the pioneering work by Aharonov 
and Bohm (AB) [2], where they proposed the first example 
of generation of topological phase acquired by an electron 
when it travels through a magnetic field-free region. This 
phenomenon, known as Aharonov-Bohm effect, has been 
the usual framework for studying properties of other phys¬ 
ical systems which lead to similar effects. The first subse- 

“e-mail: frankitalyl4@gmail.com 
’’e-mail: edilbertoos@pq.cnpq.br 
‘’e-mail: lrb.castro@ufma.br 
‘’e-mail: cleversonfilgueiras@yahoo.com.br 
‘’e-mail: diegocogollo@df.ufcg.edu.br 


quent effect was the work by Aharonov-Casher (AC) [3], 
where they predicted that the wave function of a neutral par¬ 
ticle with a magnetic dipole moment acquires a topological 
phase when traveling in a closed path which encircles an 
infinitely long filament carrying an uniform charge density. 
Subsequently, several other AB-like effects were discovered 
over the last three decades (see Refs. [4-9] ). 

An important question that we address here are the ef¬ 
fects of spin on the dynamics of topological effects. The 
first work in this context was proposed by Hagen to study 
the scattering of relativistic spin-1/2 particles in an AB po¬ 
tential [10]. He showed that, by reformulating the problem 
with a source of finite radius which is then allowed to go 
to zero, it is established that the delta function alone that 
causes solutions that are singular at the origin. He also con¬ 
cluded that the modifications in the amplitude which arise 
from the inclusion of spin are seen to modify the cross sec¬ 
tion for the case of polarized beams. Hagen has also shown 
that there is an exact equivalence between the AB effect for 
spin-1/2 particles and the AC effect [11]. This fact estab¬ 
lishes the dynamics of the AC problem. However, a pecu¬ 
liarity that Hagen has not addressed clearly in their work is 
how to find the bound states energy levels. By modeling the 
problem by boundary conditions at the origin, an expression 
for the bound state energy for the AC problem was derived 
in Ref. [12]. The method used to find these energies was es¬ 
tablished in Refs. [13, 14], and it is based on the self-adjoint 
extension method of operators in quantum mechanics. 

In the AC problem, the electric field is the one gener¬ 
ated by an infinitely long, infinitesimally thin line of charge 
along the z-axis with a charge density Xi distributed uni¬ 
formly about it, namely 

Ei=2Ai^, V.Ei=2Ai-®, (1) 

P P 
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As pointed out in Ref. [11], to ensure the exact equivalence 
between the AB and AC effects, we can not neglect the V • Ei 
term. The physical implications of this term on the dynamics 
of the particle has been quite studied in recent years [15-18]. 

Another configuration field of interest is 

E 2 = ^P, V.E 2 = A 2 , ^=0, VxE 2 = 0. (2) 

This special configuration was proposed by Ericsson and 
Sjoqvist to study an atomic analog of the Landau quantiza¬ 
tion based on the AC effect [19]. They demonstrated that the 
existence of a certain field-dipole configuration in which an 
atomic analog of the standard Landau effect occurs opens up 
the possibility for an atomic realization of the quantum Hall 
effect using electric fields. This same configuration was used 
to study the Landau levels in the nonrelativistic dynamics 
of a neutral particle which possesses a permanent magnetic 
dipole moment interacting with an external electric field in 
the curved spacetime background with either the presence 
or the absence of a torsion field [20] (see also Ref. [21]). 

In this article, we analyze the planar motion of a neutral 
particle of spin- 1/2 interacting with both electric fields of 
Eqs. (1) and (2), i.e.. 


2 Equation of motion 


We start with the Dirac equation with nonminimal coupling 
{h = c= 1 ) 




V/ = 0, 


(4) 


where ji is the magnetic dipole moment of the particle, 
is the electromagnetic tensor whose components are given 
by 

Foi = E‘ = -Ei, (5) 

— ^ijk^kt ( 6 ) 

and 




(7) 


where is the spin vector, are the components of the oper¬ 
ator 


= ( 8 ) 

As the particle interacts only with electric fields, we con¬ 
sider only (5). In the above representation, Eq. (4) can be 
written as 


E — El -f E2. 


(3) [j3 Ypi + I5m- ipj5 fEi] w = S’xj/. 


(9) 


Although the field configuration (1) has been studied in dif¬ 
ferent contexts in the literature, in our approach, we solve 
the first order Dirac equation and derive their solutions giv¬ 
ing a focus to the effects due to the spin. This analysis, in 
particular to the electric field of Eq. (1), which is responsi¬ 
ble for the AC problem, it is presented and discussed in de¬ 
tail here for the first time. We also address the second-order 
Dirac equation, which results in the Pauli equation. We make 
use of the self-adjoint extension method [ 22 ] and model the 
Hamiltonian by boundary conditions [23]. We also deter¬ 
mine an expression for the energy levels of the particle and 
compare it with the known results in the literature. 

This paper is organized as follows. In Sec. 2, we con¬ 
sider the Dirac equation with nonminimal coupling and con¬ 
struct the set of hrst order differential equations. In Sec. 3, 
we solve the first order differential equations and obtain the 
bound state solutions of the particle. The existence of these 
solutions means that the system admits isolated solutions. 
In Sec. 4, we derive the second-order equation (Pauli equa¬ 
tion) and solve it by assuming that the particle describes a 
circular path of constant radius. In Sec. 5, we analyze the dy¬ 
namics of the system with the inclusion of the r = 0 region. 
We use the self-adjoint extension method to fix the physics 
of the problem in the r = 0 region. Expressions for the wave 
functions and energies are obtained, without any arbitrary 
parameter which arises from the self-adjoint extension ap¬ 
proach. In Sec. 6 , the results of Sec. 5 are examined in the 
nonrelativistic limit. In Sec. 7, we present our conclusions. 


Following Ref. [10], we write / as 

/ = we,7)3 7^', (10) 

so that Eq. (9) becomes 

[j3 fpi + I5m- psP YEi] xj/ = i\\f, 

where £, = SijEj, Etj = —Eji- Equation (9) can be written as 

[j37-(p-)ts£)-|-j3m]'P = ^t/. (11) 

The 7 matrices are conveniently defined in terms of the Pauli 
matrices as [ 10 ] 

]37i = (7i, Pr2 = sa2, [3 = (73, (12) 

where s is twice the spin value, with s = 4-1 for spin “up” 
and s = —1 for spin “down”. Thus, Eq. (11) can be written 
as 


[ai{pi-psEi)+sa2{p2-psE2) + (y3M] (13) 

By noting that E\ = E 2 e E 2 = —Ei, as usual, we write Eq. 
(13) in polar coordinates (p, (p) 


1 ^ -frii 

^ dp *'*p id(p * \ p 
= e+“'^ (S‘-M)\i/i, 



(14) 


d 1 d 

—i^ - \-is—— -hi 

dp p id(p 


Hi 

P 


mp 




= e-“^ (<f+M)v72, 


(15) 
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where rji = 2/rAi, 772 = /rA 2 / 2 . Using the decomposition 


¥ = 



L,nfn,{p)e‘”'^ 

_ 

\¥2j 



(16) 


where m = 0 , ± 1 , ± 2 , ±3,is the angular momentum quan¬ 
tum number, Eqs. (14)-(15) provide two coupled first-order 
radial equations 


d sm + rji + 1 

+-- 

dp p 

d sm + ri] 

-b- - + 

dp p 


mp gm{p) = {S‘-M)f„{p), (17) 
mp fmip) = {S’+M)g,nip). (18) 


The factor i on the lower spinor component in Eq. (16) was 
inserted to ensure that the radial part of the spinors is man¬ 
ifestly real. An isolated solution for the problem can be ob¬ 
tained considering the particle at rest, i.e., ^ = ±M. Such 
solution for the Dirac equation in (1 -b 1) dimensions was 
investigated in Ref. [24] (see also Refs. [25, 26]). 


3 Isolated solutions and the Aharonov-Casher problem 

In order to obtain isolated solutions, let us look for bound 
state solutions subjected to the normalization condition 


^ {\fm{p)\^ + \gm{p)\^)pdp = l, 


and consider the conditions S = ±M stated above. 


3.1 Case = M 


(19) 


In this case, Eqs. (17)-(18) are written as 

d sm + 77i -b 1 

Tt—'- 

dp p 

d sm -b Til 

-—+- 

dp p 


+ ri2p 

mp 


gm{p) — 0, 
fm{p) — 2Mg^[p). 


The solutions to gm{p) and fm{p) are 

fn,{p) = 

X [ci-bC2M(T72)"“+'’‘ r (-im-77i, 772P^)] 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


where ci and C 2 are constants, and F (^—sm — r}i,r} 2 P^) is 
the upper incomplete Gamma function, obtained through the 
relation [27] 


roo 

"{a,x)= / F^^e^'dt, 91(a) > 0. 

Jx 


As rii 2 ^0, then gm{p) converges as p 0. Moreover, 
since F i^—sm — rji, J 72 P^) always diverges, then fm{p) will 
only converge if C 2 = 0 e 772 < 0. As a result, we have 

rs = ±i, 

Pi ^ 0 , 

P 2 < 0 , 

,C2 = 0, 


fn{p) 

ginip) 




(25) 


3.2 Case ^ = -M 

In this case, Eqs. (17)-(18) become 

d sm -b Tji + 1 

-7—I- 

dp p 

d sm + rji 


mp 

mp 


gm{p) — 2M/m(p), 

ftn{p) = 0. 


dp p 
The solutions of these equation are 

/,„(p) = cip"“+’)> 

g„(p)=e-2n2p2p-{“=+’)i+i) 

X C 2 -b ClM (- 772 ) -Om+’Jl + l) F [ms + 771 -b 1, -P2P^] 


(26) 

(27) 

(28) 

(29) 


As for the case S' =M, looking for solutions ( 28) and (29), 
we can see that the only square integrable solutions are 


fm (P ) 
gm{p) 


= C2 




-(snt-hTJi + 1 ) 



(30) 


In summary, note that the above results, Eqs. (25) and (30) 
are bound-state solutions of square-integrable because the 

±7)2 2 

function with 772 ^ 0 predominates over the poly¬ 

nomials p')i+"“ and p-(™+’li+P for S = M and S = —M, 
respectively. We can conclude that the presence of the A 2 is 
necessary for the existence of bound states. 

4 The quadratic equation 

The equation of second order derivative of Eq. (11) is found 
to be 

[p + s(p xE)]^i/ + P(J3(V •E)v/= (S'^-Wp-) V/. (31) 


Using (3), Eq. (31) can be written more explicitly as 


5p2 


Id \ 


p dp p2 d(p^ 


2s 


Hi 

p 


mp 


1 d 

p id(p 


(^±P 2 p) 

¥+m 

U'’<">+ 2 J 

VP / J 


L p J 


¥ 


(24) = 


= IS^-M'^) xj/, 


(32) 
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In this stage, it is worthwhile to mention that Eq. (32) is 
the correct quadratic form of the Dirac equation with non- 
minimal coupling, because the singular term from V • Ei is 
considered. 

Using the decomposition (16), the equation for fm{p) 
can be obtained 


hfm{p) = 


(33) 


with E = —M^ — 2r\2 (Pi +ms+ 1), where 


h — ho + rj I 



(34) 


is the Hamiltonian system with the magnetic moment of the 
particle pointing the positive direction of z axis, and 


ho 


dp^ 


LA. 

p dp 


,2„2 


+ ^ + ’72P • 


is Hamiltonian without the 5 function, and 


(35) 



a 


Fig. 1 Plots of the energy as a function of radius for 5 = 1 and different 
values of in. 


4.1 Particle in a ring of constant radius 

Before solving Eq. (32), an interesting case which can be 
considered here is when we assume the particle describing a 
circular motion of radius p = a = const.. In this case, from 
(32) we get 


V = m + srji. (36) 

The equation for gm{p) is obtained in an immediate way. It 
is given by 

hmip)=Eg„,{p), (37) 


1 ^ \ 
-T3-+i —+ J?ia 
a locp \ a J 


n 2 


y/= Xj/. 


(41) 


According to Eq. (16), the wave function for a particle ex¬ 
ecuting a circular motion of constant radius can be written 
as 


with E = - 2772 {in + i) + T]! — 1], where 


W = 


Qm 


eimij) ^ 


(42) 


h = 


ho-Pi 



(38) 


ho 


d^ 

dp^ 


1 d 
pdp 


22 

■-^ + ^2pL 


(39) 


and 


v = OT + i + i77i. (40) 

Equation (32) governs the system dynamics. In this dynamic, 
we must consider regular and irregular solutions, since ir¬ 
regular solutions are also physical solutions for the system 
under consideration. In other words, since we consider the 
effects of the spin of the particle and because of the field 
configuration (1), the Hamiltonian will contain a singular 
potential. We will return to this problem in Section 5. 


If xf/ is an eigenvector of O’, with eigenvalue 5 = ±1, the 
energy levels are given by 




+ s(^^ + P2a^ , (5 —+ 1 ) 


m 
a 

m-\-s 


• +5 


a 




(5 = -i) 


(43) 

(44) 


The profiles of the energy is shown in Eig. 1 for g = \ and 
some values of m. Eigure 1 clearly shows that both particle 
and antiparticle energy levels are members of the spectrum. 
Note that for positive (negative)-energy we find that the low¬ 
est quantum number corresponds to the lowest (highest) en¬ 
ergies, so that it is plausible to identify them with particle 
(antiparticle) energy levels. Also, it is noticeable that the 
Dirac energies are symmetrical about ^ = 0 and since the 
positive and negative energies never intercept one can see 
that there is no channel for spontaneous particle-antiparticle 
creation. If 772 = 0, we obtain the energy levels of a neu- 
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tral particle with magnetic moment /i in a circular path of 
constant radius. 


-M^ = \ {m + sr]i)^, (45) 

-M^ = X: (m + 5 + 577i)^. (46) 


studying the asymptotic limits of Eq. (50), we find the solu¬ 
tion 

/m(p)=P*^e-§E(p), (51) 

where (±) refers to the regular (irregular) solution, respec¬ 
tively. Substituting Eq. (51) into Eq. (50), we find 


These energies correspond to spectrum for the usual Aharonov- 
Casher effect. 


.d^F 
' dp^ 




dF{p) 

dp 




F{p)=0. 

(52) 


5 Self-adjoint extension analysis and the dynamic 
inclnding the r = 0 region 


Equation (52) is of the confluent hypergeometric equation 
type 


In this section, we solve Eqs. (33) and (37), including the 
term 5 (p) /p. In order to deal with this kind of point inter¬ 
action potential, we consider the self-adjoint extension ap¬ 
proach [28, 29]. In quantum mechanics, observables corre¬ 
spond to self-adjoint operators. However, in some physical 
systems, we deal with differential operators for which the 
Hamiltonian is not necessarily symmetrical in some region 
of the space. In such cases, the Hamiltonian is not essen¬ 
tially self-adjoint and one attempts to And self-adjoint exten¬ 
sions of the Hamiltonian corresponding to different types of 
boundary conditions. Such self-adjoint extensions are based 
in boundary conditions at the origin and conditions at infin¬ 
ity [30-32]. Erom the theory of symmetric operators, it is 
a well-known fact that the symmetric radial operator Iiq (as 
in Eq. (34)) is essentially self-adjoint for |v| > 1, while for 
IVI < 1 it admits an one-parameter family of self-adjoint ex¬ 
tensions [22], /2o,A„> where Xm is the self-adjoint extension 
parameter. Here, we will use the approach of Ref. [28, 29], 
which is based in a boundary conditions at the origin. Thus, 
all the self-adjoint extensions Iiq of !iq are parametrized 
by the boundary condition at the origin 

Ko=X„,Ki, (47) 


with 


= lim pl''l/m(p), 

p^ 0 + 


= lim —TT 

p^o+pl''l 


fm(p)-xo- 


]v| 


(48) 

(49) 


Eor Xm = 0, we have the free Hamiltonian, i.e., without the 
5 function, with regular wave functions at the origin; for 
X,„ 7 ^ 0, the boundary condition in Eq. (47) permit an 
singularity in the wave functions at the origin. Thus, by mak¬ 
ing a variable change, p = ri 2 P^^ Eq. (33) reads 




AV 

4 4772/. 


fm{p) — 0 . 


(50) 


As mentioned above, the boundary condition (47) allows us 
to look for regular and irregular solutions for Eq. (50). By 


zF''{z) + {b-z)F'{z)-aF{z)=Q. (53) 


In this manner, the general solution for Eq. (50) is given by 


fm{p)=amP^e ^F 


. - bl p 
Fb,np-^e-^F 




In Eq. (54), F{a,b,z) is the confluent hypergeometric func¬ 
tion of the first kind [27] and a„j and b,n are, respectively, the 
coefficients of the regular and irregular solutions. 

Now, we remark that Eq. (52) is equivalent to Eq. (38) 
of Ref. [33], in which the procedure for obtaining the en¬ 
ergy levels for different values of the self-adjoint extension 
parameter is given in detail. In this procedure, we use the 
boundary condition (47) together with the normalizability 
condition to obtain a relation that allows us to eliminate 
and bm of Eq. (54). Then, following Ref. [33], such condi¬ 
tion is found to be 


1 r(i + |v|) 

r(bP- 4 ) 


Equation (55) gives a contribution of the irregular solution 
to the problem. This feature comes from the fact that the 
operator Ho is not self-adjoint for |v| < 1. 

We now analyze the following points in Eq. (55): 

(i) Eor Xm = 0, case in which the 5 function is absent, 
only the regular solution contributes for the bound state wave 
function. 

(ii) Eor Xm = only the irregular solution contributes 
for the bound state wave function. 

Thus, for all the other values of the self-adjoint extension 
parameter, both regular and irregular solutions contributes 
for the bound state wave function. Analyzing the poles of 
the Gamma function in Eq. (55) together with the criteria (i) 
and (ii), we get 

= -n, for X,n = 0, (56) 

2 4772 

^ - A = -n, for Xm = (57) 

2 4772 
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Fig. 2 Plots of the energy as a function of Ai for i = 1 and different Fig. 3 Plots of the energy as a function of Ai for i = — 1 and different 
values of n and m. values of n and m. 


with n a nonnegative integer, n = 0,1,2,.... By solving Eqs. 
(56) and (57) for we obtain, respectively, for the 

regular and irregular solutions 


27)2 

2J72 


= ( 2 n + 1 + |m + sT 7 i|) + 771 + ms+ 1 , 


= (2n + 1 — \m + i 77 i I) + 771 — mi + 1 . 


Now, we consider the solution of Eq. (37). By perform¬ 
ing the same steps to achieve Eqs. (58)-(59), we obtain 




(58) 

2772 

(59) 

<^2-77/2 

2772 


= (2n -f 1 -I- |m -P i -f i77i I ) -f i (m -f i) -t- 771 — 1 , 

(62) 

= {2n -P 1 — |m-|-s-|-sr7i|) -l-s(m-l-s) -(- 771 — 1, 


As an illustration, the profiles of energy as a function of Xi 
and with spin projection parameter i = 1 and i = — 1 are 
shown in figures 2 and 3, respectively. Once again, we note 
that both particle and antiparticle energy levels are members 
of the spectrum. Also, it is noticeable that in both figures the 
Dirac energies are symmetrical about <? = 0 and, since the 
positive and negative energies never intercept, we can see 
that there is no channel for spontaneous particle-antiparticle 
creation. In this case, from the requirement of real energies 
(from equation (58)) we obtain a constraint on the minimum 
value of Xi. The parameter Xi has to satisfies the following 
inequation: 


\m + 2iisXi\ +2jj.Xi ^ 


m 2 

—;—P 2« -P 2 -P sm 

pX2 


(60) 


(63) 


with unnormalized bound state wave functions given by 

gm{p) = p^e^^E(-«,l-P|v|,p) 

+ p-^e-§E(-n,l-|v|,p), (64) 

where 


l±|v| 


E 

4772 ' 


If 771 0 in Eqs. (58)-(59) and (62)-(63), we obtain 


2772 

2772 


= {2n -p I -P |m|) -pms -P 1, 


= (2n -P I — |m|) -P mi -P 1, 


(65) 


( 66 ) 

(67) 


The unnormalized bound state wave functions are given 
by 

/m(p) = p^e^^E(-«,l-pv,p) 

+ p^^e^^F{-n,l-v,p). (61) 


and 

# 2 -m2 
2772 

#2-7^2 

2772 


(2n -pl-p|m-ps|)-psm—1, 
(2n -P 1 — |m -P i|) -P im — 1 . 


( 68 ) 

(69) 


Note that when |v| > 1 or equivalently when the 5 interac¬ 
tion is absent, only the regular solution contributes for the 
bound state wave function {bm = 0), and the energy is given 
by Eq. (58). 


These energy levels correspond to the analogue Landau quan¬ 
tization for relativistic quantum dynamics of neutral fermions 
of spin-1/2 with magnetic moment p in the field configura¬ 
tion of Eq. (2). 
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6 Nonrelativisitic limit 

Let us now examine the nonrelativistic limit of Eq. (31) by 
setting § = M + e, with e <^M. The equation to be solved 


^[p + i(MxE)]V+^M(^3(V-E)v/ = ev/. 


(70) 


Performing the same steps as for the relativistic case, we find 
the energy levels 


e = 


to [(2« + 1 + |m + i77i|) + rii +ms+ 1], 
CO [{2n + 1 — |m + i77i I) + rii+ms+l]. 


(71) 


where we have defined the frequency (O = 7)2/M. Similarly, 
we can also find the eigenvalues of Eq. (37). The result is 
given by 


e = 


CO [{2n + 1 + |m + s + sT]i|) + s(m + s) + T]i — 1], 
CO [{2n + 1 — |m + s + s77i|) + s(m + s) + T]i — 1]. 


(72) 

If T 7 i —0, we obtain the energy levels corresponding to a 
neutral fermion of spin-1/2 with magnetic moment in the 
nonrelativistic regime 


e = 


and 


e = 


CO [{2n + 1-1- \m\) — ms + 1], 
CO [{2n + 1 — \m\) —ms+l]. 


CO [{2n + 1 + |m + i|) +i (m + i) — 1], 
CO [{2n + 1 — |m + i|) +i (m + i) — 1]. 


(73) 


(74) 


Equations (73) and (74) can be compared, for example, with 
Eq. (25) of Ref. [34], in the absence of the spin element s. 


7 Conclusions 

We have solved the quantum dynamics of a neutral fermion 
with a magnetic moment ji in the presence of external elec¬ 
tric fields. We shown that the set of first order differential 
equations admit isolated solutions (E = M, and E = —M). 
This result implies new solutions for the AC problem. We 
derive the second-order Dirac equation to study the motion 
of the particle in two situations. Eirst, we assume that the 
particle describes a circular path of constant radius, and then 
analyze the dynamics in the full space, including the r = 0 
region. The inclusion of the r = 0 region states that we must 
consider the term V • Ei in Eq. (31), i.e., a singular term. We 
consider the self-adjoint extension method and show that 
the term V Ei, which results in a 5 function, has physi¬ 
cal implications on the dynamics of the particle. In other 
words, we have verified that this term contributes for the 
bound state wave function and energy spectrum, and with 
an explicit dependence on the spin projection parameter s. 


Eor Xm = 0, case in which the 5 function is absent, only the 
regular solution contributes for the bound state wave func¬ 
tion. Eor two particular values for the self-adjoint extension 
parameter, X,„ = 0 (regular solution) and A,„ = °° (irregular 
solution), the energies are given explicitly in Eqs. (58)-(59) 
and (62)-(63). In the limit Tji —0, the corresponding energy 
levels are analogous to Landau levels. In this limit, the de¬ 
pendence on s parameter is still maintained. We also have 
obtained these results in nonrelativistic limit. 
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